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Abstract 

Ward identities for Cooper pairs are derived. These give consistent 
description of electronic curent vertex and thermal current vertex. 
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The quest for the correct expression of the thermal current vertex for Cooper 
pairs has a long history but the present status is still controversial. 

A phenomenological expression on the basis of Ginzburg-Landau (GL) the- 
ory jTl[2] is naively expected to be reliable. The GL theory relates the electronic 
current vertex J e and the thermal current vertex for Cooper pairs as 

in the limits of long wavelength and low frequency where incoming and outgoing 
Cooper pairs have bosonic thermal frequencies iui m and iui m + iw v and e is the 
charge of an electron (e < 0). 

On the other hand, the most recent works [31 0] in this field relate these two 
vertices as 

j Q = l^u m + ^r. (2) 

However the origin of the factor 2 has not been explained convincingly. 

In this Short Note we try to obtain the correct relation between J e and J® 
on the basis of Ward identities. 



•APPENDIX has been added in v2. 
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First we review the derivation [5] of the Ward identity for electronic current 
vertex. The vertex function is defined as 

A« y, z) = (T T {j*(z)^(x)i;\(y)}), (3) 

and under the charge-current conservation Q3u=o af - -^ 2 ) = ^) ^ s divergence 
is transformed into 

3 a 

- l Y,Q^ A >>y> z ) = {Tr{\jl{z)^{x)]^{y)})8{T z -T x ) 

where (■ • •) represents the thermal average, T T is the time-ordering operator 
with respect to the imaginary time r z = izo and z = (z, Zq) with the coordinate 
vector z = (zi,Z2,zz) and the real time zq. Here ipf(x) and ip^(y) are annihi- 
lation and creation operators of f-spin electron. The zeroth component of the 
electron current jo(z) is given by Jq(z) = eij}t(z)ipf(z) + eip\(z)i/)i(z). Using the 

commutation relation [io(z), = e ^\{ z ) an d \jo(z), ipt( z )} = — e- 0t( z ) an d 

introducing the Fourier transform, we obtain the Ward identity 

3 

]T fc M r« (p, fc) = eG- 1 ^) - eG^fr + fc), (5) 

for the electronic current vertex where A e ^(p,k) = G(p)T e (jp, k)G(p + k) and 
G(p) is the electron propagator with four-momentum p — (p,po). The zeroth 
components po and ko are fcrmionic (po = ~i£n) and bosonic (fco = — iu> u ) 
frequencies. 

The extension of this Ward identity to the case of Cooper pairs is straight- 
forward. In the following we consider the Cooper pair of s-wave pairing in 
the case of local attractive interaction. Since we discuss the normal metallic 
phase (T > T c ), the propagator of Cooper pairs is a fluctuation propagator. 
Replacing ip^\-(x) by ?pi(x)ipf(x) and ip\(y) by ip^(y)ip^(y) and using the com- 
mutation relation \j$(z), %l)\{z)^\{z)] = 2e^(z)^\{z) and [j^(z),^ i (z)ip- t (z)} = 
—2eipi(z)ip^(z), we obtain the Ward identity 

3 

k^l(q, k) = 2eD-\q) - 2eD-\q + k), (6) 

for Cooper pairs where is the counterpart of and D{q) is the Cooper-pair 
fluctuation propagator with four- momentum q = (q,qo) whose zeroth compo- 
nent qo is a bosonic frequency (qo = —iuj m ). It should be noted that the factor 
2e represents the charge carried by a Cooper pair and is automatically taken 
into account by the commutation relations. This point was missed in an early 
guess [6] of the Ward identity for Cooper pairs. 
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Second we review the derivation [7J of the Ward identity for thermal current 
vertex. Since the essence of the derivation is the same as the electronic current 
vertex, this review is short. In the limit of vanishing external momentum, k — > 0, 
the Ward identity is obtained as 

3 

]T k^(p + k,p)= PoG-^p + k)- (po + fco)G- 1 (p), (7) 
n=o 

for electrons. Here po and pa + fco result from the Fourier transform of the time- 
derivative of annihilation and creation operators of t-spin electron. The time- 
derivative results from the commutation relation between the zeroth component 
of the thermal current and annihilation or creation operator, since the zeroth 
component for k —> is the Hamiltonian of the system. 

The extension of this Ward identity to the case of Cooper pairs is also 
straightforward. The Ward identity for thermal current vertex for Cooper pairs 
is 

3 

J2 k^(q + k,q) = q D-\q + k) - (q + k )D- l (q). (8) 

fj,=Q 

Finally our Ward identities eqs. (6) and (8) are consistent with the GL result 
eq. (1). 

This work arose from discussions with Kazumasa Miyake, Yukinobu Fuji- 
moto and Shinji Watanabe at Osaka University. 

APPENDIX 

In this APPENDIX the derivation of the Ward identity for Cooper pairs is 
explained in detail. 

1 Introduction 

The Ward identity for electric current vertex is explained in Schrieffer's text- 
book [5]. 

The Ward identity for heat current vertex is derived in Ono's paper [7J. 

If we understand these discussions well, we can reach the Ward identity for 
Cooper pairs with little effort. Thus I review these works first and cast the 
results into those for Cooper pairs. 
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2 Ward Identity for Electric Current Vertex 

First we review the derivation [5 of the Ward identity for electric current vertex. 
We mainly discuss the case of zero temperature and cast the zero-temperature 
result into that of finite temperature. 

We consider the three-point function (fi= 1, 2, 3, 0) defined as 

A*(x,y,z) = (T{jZ(z)^(x)ip\(y)}), (9) 

where (• ■ ■} represents the expectation value in the ground state, T is the time- 
ordering operator and z = (z,zq) with coordinate vector z = (zi,Z2,Z3) and 
time zq. Here ipt( x ) an d V'-j-(y) are annihilation and creation operators off-spin 
electron. The electric current jfj obeys the charge-conservation law 

£w/>m=°- (io) 

Especially the electric charge jg is given by 

f (z) = e^\{z)^(z) + e^l(z)^(z), (11) 

where e is the charge of an electron (e < 0). The time ordering of three operators 
results in the summation of 3! terms as 



^(x,y,z) = (j e ^(z)i; t (x)ip\(y))0(zo 


- x o )0(x o 


- yo) 


- (jZ(z)4(y)i> t (x))6(z - 


- VoWiVo - 


- Xq) 


+ (^(aOj«(*W{(y)>0(s o 


- z )9(z Q - 


- yo) 


- (ip\(y)j;(z)i> t (x))e( yo - 


- z )6(z - 


- x ) 


+ (i> t (x)ip\(y)j;(z))6(x 


- yo)d(yo ■ 


- z ) 




- x )9{x 


- zo), 



where 0(x) is the unit step function. Thus the time-derivative of A^ results in 

— A-l(x,y,z) = S(z -x )[9(x - y )([jo(z), tfo(x)]ip\(y)) 

~9(y a -x a )(4>\(yM(z),^(x)})) 
+S(z - yo)[9(x - j/oK^tMboC*)) ip\(y)}) 
-6{y Q -x )([f Q {z)^\{y)]^{x))) 

HT{?§^M*)4(y)})- (is) 
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Using again the time ordering the divergence of is expressed as 
3 1 

J2 jrK^ V> z ) =(T{{f (z),^(x)H(y)})6(zo - x ) 

fj,=0 ^ 

+(T{^(xM(z),Tp\(y)]}}5(z -y ) 

+( T |E^^^t(^K(y)|). (14) 

The last term on the right-hand side vanishes due to the charge-conservation 
law, ea. (|10p . Only equal space-time commutation relations are non- vanishing, 

{f (z)^\(z)}=e^\(z), (15) 

and 

[j e (z),^ t (z)] = ^ t ( 2 ), (16) 
so that the non- vanishing contribution becomes 
3 a 

J2 y, z) = - e(T{M^\(y)})S 4 ^ x) 



/i=0 



+ e{T{^{x)^\{y)})5\z-y). (17) 
Introducing the electron propagator G(x, y) as 

G(x,y) = -i(T{^ f (x)^\(y)}), (18) 
this relation is written into 

3 q 

V — A e (x, y, z) = -ieG(x, y)S 4 (z - x) + ieG(x, y)5 i (z - y), (19) 

OZn * 

Assuming the translational invariance we set y — and introduce the Fourier 
transform as 

A£(p, k) = J d 4 xe~^ J d 4 ze- lkz (T{j;(z)Mx)^\(0)}), (20) 

where the four-momentum is defined as p — (p,po) and k — (k,ko). The left- 
hand side of eq. P^|) is evaluated as 

t^>M = j ^ iPX j ^^pk^k), (21) 
and the right-hand side is transformed as 

^-»./^--(- G M)^- l) + G (M)^)) 

= -G(p + k)+G(p), (22) 
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where 



G(p) = J d 4 xe- ipx G(x,0). (23) 



Therefore we obtain 

3 

Y / k l ,Al(p,k) = eG(p)-eG(p + k). (24) 

fl=0 

The vertex function rf. is introduced as 

A£(p, AO = iG(p) • r£(p, *)• iG(p + *), (25) 

in accordance with the definiton of the Green function, eq. (|18l) . Then the Ward 
identity for the electric current vertex is given by 

3 

Y, W(P. fc ) = eG ^iP) ~ eG ^(P + k )- ( 26 ) 

Since the Fourier transform is introduced as 

px =pixi +p 2 x 2 +P3X 3 - et, (27) 

where e is the energy and t is the time, xo — t and po = —e (and in the 
same manner ko = —lu with uj being the energy of the external field) in the 
zero-temperatute formalism. 

Here we check the limiting case of eq. (|26|) . If we replace the full Green 
function G{p) by the free Green function Go(p) and set ko = 0, we obtain 

£W(P'*) = -^ (28) 

where the free dispersion eg = p 2 /2m is employed with m being the mass of 
electron. This relation means 

rj(p,0)=e«„, (29) 

where the right-hand side is the proper electric current vertex with the electron 
velocity v = p/m. 

In the finite-temperature formalism we employ the time-ordering operator 
T T and consider the three-point function 

A;(x, y, z) = (T T {j;(z)^{x)^\(y)}), (30) 

where the real time zq and the imaginary time r z is related by r z = izq and (• • •) 
represents the thermal average. Taking the charge-conservation law, eq. (|10[) . 
into account we obtain 

~ l E nrK^ V' z ) =( T AUo(^ Mx)]4(y)})s(r z ~ r x ) 

l-L — ^ 



+ {Tr{^{xM{z)^\{y)]})8{T z - t v ), (31) 



G 



instead of eq. ((HJ). Here we have used the relation 

d . d 

dr z dz ' 



(32) 



The finite-temperature vertex function T e ^ is introduced as 

Klip, k) = [~G{p)\ ■ r£(p, k) ■ [-G(p + k)}, (33) 

in accordance with the definiton of the thermal Green function 

G(x,y) = -{T r {^(x)^\(y)}). (34) 

Then the resulting Ward identity is the same form as eq. (j2"6")l in the case of zero 
temperature. For the finite-temperature Ward identity the zeroth component 
of the four-momentum is po = —is n with fermionic thermal frequency e n and 
fco = — iu) v with bosonic thermal frequency w„. 

The extension of this Ward identity to the case of Cooper pairs is straight- 
forward. In the following we consider the Cooper pair of s-wave pairing in the 
case of local attractive interaction. Replacing 4"[{ x ) by ^(x) — ifji(x)ipf(x) and 
ijj^(y) by ^(y) = ip^.(y)ip^(y) we consider the three-point function as 

Mfa,y,z) = (T{jfc)*(x)&(y)}}, (35) 

where ^(x) and &'(y) are annihilation and creation operators of a Cooper pair 
which has a bosonic character. Using the commutation relation 

[f Q (z),^(z)}=2e^(z), (36) 

and 

[f (z)^(z)} = -2e^(z), (37) 
the divergence of is expressed as 

-£-M;(x,y,z) = - 2e(T{*(x)&(y)})S\z - x) 

+ 2e(T{*(x)^(y)})6\z-y), (38) 



by repeating the same calculations as eqs. (|12|) . (|13|) and ([14]) . Here the differ- 
ence between fermion and boson is handled solely by the time-ordering operator 
T so that the expression of the divergence is common to fermion and boson. In- 
troducing the Cooper-pair propagator D{x,y) as 

D(x,y) = -i(T{9(x)9\y)}), (39) 
we obtain the Ward identity 

3 

^ k^(q, k) = 2eD- 1 {q) - 2eD-\q + k), (40) 
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for Cooper pairs where is the counterpart of and D(q) is the Fourier 
transform of D(x,0) with four-momentum q — (q,qo)- 

Although the above derivation for Cooper pairs is formulated at zero temper- 
ature, eq. (|40p also holds at finite temperature with q$ being a bosonic thermal 
frequency (go = —i&m)- We are mainly interested in the normal metallic phase 
(T > T c ), the Cooper-pair propagator is a fluctuation propagator in this case. 

It should be noted that the factor 2e represents the charge carried by a 
Cooper pair and is automatically taken into account by the commutation rela- 
tion. 

3 Ward Identity for Heat Current Vertex 

First we review the derivation [7] of the Ward identity for heat current vertex. 
We consider the three-point function AJj defined as 

A?(x,y,z) = (T{jQ(z)^(x)i;\(y)}}, (41) 

where is the heat current. The heat current obeys the energy-conservation 
law 

Assuming the translational invariance we set 

A^(x,y,z) = A^(x-y,z-x), (43) 
so that the Fourier transform becomes 

d 4 ze~ ikz J d A xe- ip ' x J d 4 ye ipy AQ{x, Vl z) = AQ{ P ,p-k)(2nf5 4 {-k-p'+p), 

(44) 

where 

A«(p,p-fc) = / d 4 (x-y)e-^-y^ [ d 4 (z~x)e- lk ^AQ(x-y,z-x), (45) 



and S(— k — p' + p) represents the conservation of four- momentum. The diver- 
gence of Ajf becomes 



3 d 



£ —AQ(x,y,z) ={T{\j$(z),MxM(y)})S(zo-xo) 

fji—0 r 

+{T{^(x)\jQ(z),i>\(y)]}}S(z - yo), (46) 
under the energy-conservation law, eq. (|42|l . The left-hand side is evaluated as 

(47) 
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so that the Fourier transform satisfies 

3 



-ik (z -x ) 



ju=0 



^2ik^(p,p-k) = J d(x - y )e-^ x <>-y°) J d(z -x )e- 

{T{\j%(xo),a p _ £ (x )]a]fao)})5(zo - x ) 
+ (r{a-_ fc -(x )[j?( 2 /o),4(yo)]}}<5(^ - Vo)), (48) 
where 

tf(z) = ]T (*o), ^r(x) - e *'*°? (*<>)> ) = E e " ip "H(^)' 

k p 1 P 

(49) 

with 

Then we evaluate the commutation relations. For such a purpose we intro- 
duce the Hamiltonian density h(z) for an isotropic system as h(z) = /i kln (z) + 
h lnt (z) at time Zo where 

^ kin ( f ) = ^ E v ^ • v ^(*)> ( 51 ) 

CT 

and 

^(3 = ^EE / (52) 

a a' J 

with a general interaction strength V(2 — z*) of two-body interaction and spin 
a, cr' =t)4- Within this paragraph the time of all operators is zq. This Hamil- 
tonian density h(z) is the zeroth component of the heat current j®{z). The 
Fourier component 

% = = j dze-^(z), (53) 

is decomposed as jS = jh ln + jl nt and 

it - i E(^ - *o • + b U b ?) ' (54) 

operators of J,-spin electron. Comparing 



where bp and 61, - are the Fourier components of annihilation and creation 



r .kin 1 _ (P~k)-P r ki „ +1 - (? ~ ' P t ,55} 
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with 

we can use the replacement 

lit^-ti^litv^ ft.«J-sl- (57) 

in the limit of vanishing external momentum, fc — > 0. The similar argument 
holds for [j?*, and [j g nt , a|r] so that we obtain 

<V-g] [#> « P i b|, 4 (58) 

since jS = H where 

J h(z)dz, (59) 



'k=0 

H 



is the Hamiltonian. 

Thus the equation of motion 



d d 
[H,ap{x )} = -i—a p {x ), [H, a^_ % {y )] = - i g^ a l^(yo)> ( 60 ) 

can be applied to cq. ([4"5)l and the result is 

Jjk li A$(p,p-k) = J d(x - yo)e-* Mx °- yo) J d(z - x ) e -^-^ 
/ d d \ 



M=0 

X 



where 

Gp(x - yo) = -i(T{ap(xo)al(y )}) ■ (62) 

Employing 

^Gfco -yo)=J ^j p '° (xo - yo) ip'oGM), (63) 

the first term on the right-hand side of eq. ([BTjl becomes ipoG(p) where G(p) = 
Gp(po)- The second term is integrated in the same manner and becomes —i(po — 
k )G(p — k). Then shifting the four-momentum we obtain 

3 

J2 fc/.A?(P + k,p) = (po + k Q )G(p + k)-p G(p). (64) 

jll=0 

This relation is converted into the Ward identity 

3 

J2 k^(p + k,p)= Po G-\p + k)- (p + fco)G- 1 (p), (65) 
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for heat current vertex. Here we have used the same relation A®(p + k,p) 
iG{p + k) ■ TJ3 (p + k,p) ■ iG{p) as eq. 



In the case of finite temperature the same relation as eq. (|65|) holds. 
If we only consider the local interaction, the derivation becomes very sim- 
ple [5] as follows. Since we are mainly interested in BCS-type local interaction, 
such a derivation is sufficient for our purpose. Taking into account the fact, 
jo(z) — h(z), eq. (|16"|) is equivalent to 

3 r, 

7^A«(x, y, z) =(T{[h(x), ^{x)]4,\{y)})8\z - x) 

fl—0 ^ 

+{T{^(x)[h(y),^\(y)]})5 4 (z-y). (66) 

If the interaction among electrons is local, then [h(x),^(x)] = [H, ip^(x)] and 
[h{y)^\{y)\ = [H^\{y)\. Thus using 

[h{x), ^{x)\ = -i-lL^x), [h{y)^\{y)] = -i-S^\{y), (67) 
eq. (1551) becomes 

V ^-A«(x, y, z) = ^-G(x, y)5\z - x) + ^-G(x, y)8\z - y), (68) 
oz M " dx dy 



which has the similar structure as eq. (|19p . Assuming the translational invari- 
ance we set G(x, y) — G(x — y) and introduce the Fourier transform 

^Gi X -y) = f^e^^G^ (69) 
as eq. (l63l) . Performing the Fourier transform in eq. (|45|) we obtain 

3 

53 k^( P ,p- k) = PQ G(p) - (p - k )G(p - k), (70) 

(ii=0 

which becomes eq. (|64[) by shifting the four-momentum. 

Here we check the limiting case of eq. (|65|) . If we replace the full Green 
function G{p) by the free Green function Go(p) and set fco = 0, we obtain 

53 * p r«(p + A,p) = ~*-(?+|), (7i) 

^ — ' ^ ml 

it=i 

and via po = — e this relation means 

rQ(p,p)=ev IM (72) 
where the right-hand side is the proper energy current vertex. 
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The extension of this Ward identity to the case of Cooper pairs is also 
straightforward. The Ward identity for heat current vertex for Cooper pairs 
is 

3 

Y, + k >l) = 1oD-\q + k)- (q + k^D-^q), (73) 

n=o 

where is the counterpart of T®. 

It should be noted that the factor q represents the energy carried by a 
Cooper pair and is automatically taken into account by the commutation rela- 
tion. 

4 Conclusion 

I have explained the derivation of the Ward identity in detail. I hope that you 
can save time in understanding it by this APPENDIX. 

Our main results of Ward identities for Cooper pairs are eqs. (|4"0")l and ([731 
whose counterparts for electrons are eqs. (|26|) and ([65]). The charge and energy 
carried by Cooper pairs are properly taken into account by the commutation 
relations. 
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